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On  the  velocity  distribution  of  turbulent  flow  behind  a 
system  cf  thin  cylindrical  rods 
By  R.  Gran  Olsson* 

1.  Introduction.  In  an  earlier  paper  the  author  has 
treated  the  problem  of  the  velocity  distribution  of  turbulent 
flow  behind  a system  of  rods  spaced  at  equal  distances  from  each 
other  [1],  For  the  shear  stress  of  turbulent  flow  the  expression 
in  terms  of  the  mixing  length  due  to  L.  Prandtl  was  assumed, 
where  the  mixing  length  at  a certain  distance  from  the  rods  may 
be  put  equal  to  a constant  t [2],  The  velocity  distribution  was 
given  by  a non-linear  differential  equation  of  the  second  order, 
the  solution  of  which  can  be  expressed  by  quadrature.  The  purpose 
of  this  note  is  to  show,  how  a solution  of  the  differential  equa- 
tion can  also  be  expressed  in  terms  of  the  elliptic  functions  of 
Weierstrass  [ 3 ] • 

2,  The  differential  equation.  We  consider  a system 
of  cylindrical  equidistant  rods  of  infinite  length  in  a plane 
perpendicular  to  the  direction  of  flow  in  e stream  of  originally 
constant  velocity.  The  problem  is  to  calculate  the  distribution 
of  velocity  in  the  region  behind  the  rods  by  means  of  the  ex- 
pression of  L.  Prandtl  for  the  shear  stress  in  turbulent  flow. 

The  mixing  length  will  also  be  constant  in  a plane  parallel  to 
the  plane  of  the  rods  because  of  their  uniform  spacing. 

The  coordinate  system  may  be  chosen  as  follows: 

The  y-axis  parallel  to  the  plane  of  the  rods  and  perpendicular 
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to  their  axes,  the  x-axis  through  the  axis  of  one  of  the  rods  and 
parallel  to  the  direction  of  the  undisturbed  flow  (Fig.  1).  The 
components  of  the  velocity  nay  be  noted  by  u and  v,  the  velocity 
of  the  undisturbed  flow  by  U.  Behind  the  rod  system  at  a suf- 
ficiently large  distance  the  components  of  velocity  are 

u = U - ulf  v = v-p  (1) 

where  u-^  and  v^  are  assumed  to  be  small  compared  with  U.  In 
establishing  the  equations  of  motion  all  terms  of  higher  order 
of  magnitude  are  dropped ; furthermore  the  pressure  term  is  neg- 
lected. This  assumption  can  be  justified  bv  means  of  the  solu- 
tion obtained. 

The  equation  of  motion  in  the  x-direction  is 


3u  , 

u aS  + v 


3u  _ 1 3a 

ay  ” p 3y  * 


(?) 


Introducing  u = U - u^,  v = v^,  where  u^  <<U,v^«  U,  we  obtain 
as  a first  approximation  the  equation 


3ui  x 9t 

9x  p 9y 


(2a) 


Because  of  the  periodicity  of  the  rod  system,  it  seems  convenient 
to  assume  the  following  expression  for  u-^: 


u^  = Af (x)cos  ay,  (3) 

where  a has  to  be  chosen  in  such  a manner  that  the  period  of 
agrees  with  the  distance  X of  the  rods,  i.e.  a=  2nA» 

The  expression  for  shear  stress  according  to  the  theory 
of  mixing  length  is 
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The  differential  equation  (2a)  wjII  then  be  transformed  into 


9u, 


au. 


a£ 


ui 


SIT  = “k  3F 

where  k = 2-t^/U  is  a constant  parameter. 


(2b) 


The  solution  of  the  differential  equation.  If  we 
make  the  substitution 

U1  = Fi(y'  (5) 

with  the  condition  F^(0)  = 1,  A may  be  interpreted  as  the  ampli- 
tude of  the  velocity  distribution,  and  we  get  the  following  dif- 
ferential equation 

F,(y)  = A kF«(y)F"(y)  (2c) 

x oil 

Assume  now  y = cq,  where  c is  a constant  of  the  dimension  length, 
we  obtain  the  following  equation 

F (cq)  = ^ F»(cn)F»(cn)  (2d) 

1 c2  1 1 

where  on  the  right  side  the  differentiations  are  made  with  respect 
to  the  variable  q . As  c is  arbitrary  we  choose  it  so  that 


A k/c2  = 1 
o 

and  obtain  the  differential  equation 

F1  = F1*F1 

which,  multiplied  by  F|,  results  in 


F,F'  - F'dF" 


1 1 


1 1 


(2e) 


(2f) 


and  once  integrated 
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(2g) 


where  C is  a constant  of  integration,  which  can  be  determined  by 


the  condition  that  for  y = 0 (rpO)  F^  = 1,  and  Fj  = 0 or  C = 1/2, 
With  this  value  of  C we  get 


Fl3  = ~ 2(1  ' Fl) 


(2h) 


From  the  substitution 


2 i o 

= kz  + 1, 
1 ’ 


dF^  - (~z?  dz/F^  = 6z2  dz/(l  + 4z^)' 


(6a) 


the  independent  variable  q is  given  by 


n = <6)4 


2/3  P°  Z dz 


3+1 


which  is  an  elliptic  integral  of  V/eierstrass  type.  To  get  a more 


convenient  expression  for  q and  F,  we  introduce  the  parameter  w 


by  means  of  the  elliptic  integral 


z 


where  the  inverse  function 


z = p(v) 


is  the  elliptic  p-function  of  V.'eierstrass,  Fy  comparison  of 


the  integrals  (7)  and  (8),  we  see  that  in  our  case  g2  = 0,  g^ 
From  (8)  and  (9)  one  obtains 


n = (6)2//3  ' 0 p(w)dw  = (6)2//^  [p  (w)  - p(w  )]  (7e) 

■J  w 0 
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where  w corresponds  to  z = 0.  u to  an  arbitrary  value  of  ?.  and 
o 

p(u)  to  the  elliptic  p -function  of  Wei erstracs  defined  by  the 


integral 


p (v  ) = - p(w)dw 


The  function  F^(n)  of  velocity  distribution  is  given  by  Eqs.  (6) 


and  (9) 


F^(w)  = {~4r>^(w)  + lj^2 


(6a) 


On  the  other  hand  the  derivative  of  the  weierstrass  p-function 


is  given  by 


p’(w)  =<j4p^(w)  - g2P(w)  - z^\> 


(10) 


so  that  the  velocity  distribution  is  given  by 


F1(w)  = p*  (w) 


(6b) 


Thus  we  have  a representation  of  the  variables  F-,  and  n by  means 
of  the  common  variable  u,  Tht  result  of  the  numerical  calculation 
is  shown  in  Fig,  2 and  is  obtained  by  using  the  tables  of  the 
Jacobian  elliptic  functions  by  t,  M.  Milne-Thomson  [4],  The  p - 
function  of  Weierstrass  can  be  expressed  by  the  Jacobian  elliptic 


function  in  the  following  manner 

l+cn(2wH1/2) 

p (w)  = e,  + H T7p7> 

2 l-cn(2wH1/£?) 


(11) 


r)'(w)  = **4h 


3/2  sn(2wH1/2)dn(2wH1/2) 
^ l-cn(2wH1/2)}2 


(11a' 


where  sn,  cn  and  dn  are  the  Jacobian  elliptic  functions, 

H2  = (eo-e, ) (e  -e. ) = 2e2  + -2— 

2 1 2 3 2 3e2 
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and  e.,  , arc  the  roots  of  the  equation 

4w3  - g2w  - = 0, 

8o  being  the  real  root,  whereas  e,  and  are  complex.  In  our 
*-  J-  3 

case  e9  = -I/\^+  = -0,6300, 

The  independent  variable  q is  expressed  by  the  integral 

of  the  p -function,  i.e, 

w 

'j  1 4-^n  t ywM — — 1 

(7b) 


,2/3  ^"°r  . l+cn(2wHl/2)  , . 

q = 6 J [ e0  + rf r7rT  J uw 

I w 2 l-cn(2wH-^2) 


The  result  is  in  close  agreement  with  the  results  obtained  by 
simple  quadrature  of  an  equation  corresponding  to  (7). 

The  velocity  distribution  in  the  direction  parallel  to 
the  plane  of  the  rods,  but  perpendicular  to  their  axes  is  given 
by  the  equation  of  continuity 


3u- 


0V 


ax  ay 

where  according  to  Eq.  (5) 


I = 0 


(12) 


~ = -2A  x‘2F,(y)  (lia) 

ax  oi 

and  therefore 

= cAox~2  ^F^cq  )dq  (12a) 

We  now  use  the  expressions  of  F-^(cq)  from  Eq.  (6b),  dq  from  Eq, 
(7)  and  z = p (w)  from  Eq.  (9),  obtaining 

v.,  = 6*"^  Acx~2  J'p(w)  p?(w)dw 

= C^)1^  Acx”2  [p^(w)  - p2(w^)  ] 

Expressing  the  p -function  by  the  Jacobian  elliptic  functions  we 
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get  a velocity  distribution  in  close  agreement  with  that  found 
by  simple  quadrature  (Fig.  3). 

4,  Summary.  The  velocity  distribution  behind  a system 
of  cylindrical  equidistant  rods  is  found  by  using  the  expression 
of  L.  Prandtl  for  the  shear  stress  in  turbulent  motion  (mixing 
length  theory).  The  solution  of  the  differential  equation  is 
given  in  terms  of  the  elliptic  functions  of  h’eierstrass , which 
for  the  numerical  calculation  are  transformed  to  the  elliptic 
functions  of  Jacobi.  The  results  are  in  close  agreement  with 
those  obtained  earlier  by  simple  quadrature. 
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